arXiv:1505.02762v2 [cond-mat.str-el] 5 Oct 2015 


Protection against Spin Gap in 2-d Insulating Antiferromagnets with a Chern-Simons 

Term 

Imam Makhfudz and Pierre Pujol 

Laboratoire de Physique Theorique-IRSAMC, CNRS and Universite de Toulouse, UPS, F-31062 Toulouse, France 

(Dated: October 6, 2015) 

We propose a novel mechanism for the protection against spin gapped states in doped antifer¬ 
romagnets. It requires the presence of a Chern-Simons term that can be generated by a coupling 
between spin and an insulator. We first demonstrate that in the presence of this term the vortex 
loop excitations of the spin sector behave as anyons with fractional statistics. To generate such term, 
the fermions should have massive Dirac spectrum coupled to the emergent spin field of the spin sec¬ 
tor. The Dirac spectrum can be realized by a planar spin configuration arising as the lowest-energy 
configuration of a square lattice antiferromagnet Hamiltonian involving a Dzyaloshinskii-Moriya in¬ 
teraction. The mass is provided by a combination of dimerization and staggered chemical potential. 

We finally show that for realistic parameters, anyonic vortex loop condensation will likely never 
occur and thus the spin gapped state is prevented. We also propose real magnetic materials for an 
experimental verification of our theory. 


I. INTRODUCTION 

The discovery of high Tc superconductivity in Cuprates 
led to a flurry of new ideas such as the idea of spin liquid 
state as one of the possible explanations for the emer¬ 
gence of superconductivity from Mott insulating par¬ 
ent compounds [1] and the importance of the physics of 
doped antiferromagnets [5]. An exotic spin liquid state 
from Heisenberg types of spin models involving break¬ 
ing of discrete symmetries such as parity or time-reversal 
symmetry, the so-called chiral spin liquid, has been pro¬ 
posed [3]. The effective low energy theory of such chi¬ 
ral spin state normally involves a topological term called 
Chern-Simons term [4]. It was originally studied in parti¬ 
cle physics [S] and mimicks the fractional quantum Hall 
effect (FQHE) [5] where such term appears as the low 
energy effective theory in the bulk [7j. 

The Chern-Simons term can be generated by a 
fermion-gauge field coupling when the fermion is inte¬ 
grated out. One obtains a fermion determinant which 
gives the Chern-Simons term as the action for gauge field 
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in Euclidean space-time upon perturbative expansion [9] , 
where Nf is the number of fermion flavors, e is the gauge 
charge, is the Dirac fermion mass, and is totally 
antisymmetric tensor. The Chern-Simons action depends 
only on the overall sign of the mass rather than its mag¬ 
nitude even though the mass must be nonzero for the 
expansion to make sense. 

In quantum magnetism, the magnetization curve may 
show the presence of plateaus. It is well understood that 
the plateau state corresponds to a state with gapped 
magnetic excitations while outside the plateau, the spin 
sector is gapless [10]. In doped antiferromagnets 
this leads to distinct natures of interaction between the 


fermions [12] . A state with gapped magnetic excitations 
is normally associated with preserved continuous sym¬ 
metry with no long range magnetic order in the direc¬ 
tion transverse to the applied magnetic field m- In this 
work, we show that in the presence of a Chern-Simons 
term one can have either one of two scenarios: i)A spin 
gapped state occurs as the analog of the fractional quan¬ 
tum Hall effect (EQHE) with its chiral edge states ii)The 
system is protected against such spin gapped state. In 
this paper, we consider 2-d antiferromagnets on square 
lattice, as we did in our previous work [la¬ 
in this paper, we also propose a way to realize a chi¬ 
ral theory as an effective low-energy theory of spin sys¬ 
tems explicitly rather than spontaneously by consider¬ 
ing doped antiferromagnets on the square lattice. The 
fermions hop on the lattice on top of a pre-existing mag¬ 
netic background. Because of a strong Hund coupling, 
the spin of the electron must be parallel to that of the lo¬ 
cal spin within the adiabatic approximation. The result¬ 
ing dynamics is well-described by effective tight-binding 
Hamiltonian [T3j 


H = —t''^^{Utr\Utr')cl,Cr + h.C. (2) 

r,r' 


where the spin sector modifies the hopping integral via 
the overlap of spin coherent states |Hr) between nearest- 
neighbor sites, a mechanism well established from the 
studies of anomalous Hall effect [T3] and doped antifer¬ 
romagnets m- This coherent state overlap has two ef¬ 
fects: providing the background flux for the fermions and 
the effective spin field that couples to the fermions. The 
spin sector should take an appropriate classical lowest- 
energy configuration which produces a staggered tt-Aux 
state [16] with dimerization and staggered chemical po¬ 
tential needed to get a massive Dirac fermion spectrum. 
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II. FIELD THEORY 

We use the semiclassical spin path integral approach 
[To] and start with the Euclidean effective action of 2-d 
doped antiferromagnet 


Euclidean space-time, we obtain the following effective 
Lagrangian for Jy 


S^ = j cPx J dT^{dr(l)f + dr(j) 

(3) 

^ - eA^) + m^]ip (4) 

taking the form of an XY model describing quantum 
fluctuations around a classical lowest-energy configura¬ 
tion specified by S = ^(sinflo cos ^Oj sindg sinc/ig, cos 6 >o) 
plus a Berry phase action whose form follows from the 
spin configuration we will describe. The (j) is the quan¬ 
tum fluctuating phase angle field around (^o |inj and is 
conjugate to the momentum operator If; [0^, Hr'] = f^rr' 
where If can be defined in terms of 6 [a- The Kr,Kr 
are stiffness coefficients of the spin sector. The ip, ip are 
Dirac spinors representing electrons that couple to the 
spin sector’s (p field via the pseudo-gauge field given as 
A^ = dfj,(p,H = T,x,y. Here two comments are in or¬ 
der. First, the definition of the vector field looks as 
if it were a pure gauge. However, one has to keep in 
mind that the field (p may have vorticity, which turns out 
to be at the origin of a non-zero flux for Second, 
the action is not gauge invariant and as such the theory 
is not a genuine gauge theory. We use nevertheless the 
terminology of pseudo-gauge field for which we also 
refer to as spin field, because of the way it couples to 
the charge degrees of freedom. When the fermions are 
integrated out, we get the Chern-Simons term due to the 
fermion-spin field coupling given in Eq. 0 . A similar 
final action but without Berry phase has been studied in 
a different context Hi!- 

Integrating out the massive Dirac fermions in Eq. 
we obtain a Chern-Simons term in terms of the phase 
field (p 

Ccs = i^t^''^d^cPd,dx<P (5) 

ZTT 

with K = e^/4 {Nf = 1). We first investigate the effect 
of this Chern-Simons term added to the action Eq.([^ 
by applying a duality transformation which introduces a 
Hubbard-Stratonovich auxiliary vector field [10] and 
re-expresses the full action as 


-|-f27r f b^{x)j!p,{x) (7) 

J X 

where = \{^)er^ivXu- As in [TU], the above result was 
obtained under the assumption that only configurations 
with closed vortex loops do contribute to the low energy 
physics. The effective magnetic field through a vortex 
loop is given by 


( 8 ) 

We note that we get an extra term /cJy to the effec¬ 
tive magnetic field coming from the Chern-Simons term. 
The important consequence of this is that the resulting 
effective Berry phase for a vortex loop now has a total 
contribution coming from ordinary Berry phase term and 
Chern-Simons term and is given by 

4^vortexloop — other vortex loop (9) 

where A is the area of a vortex loop and equals integer 
multiple of a^, q is the vorticity of the vortex loop, and 
^’other vortex loop is the flux of Other vortex loop, nonzero 
only if that other vortex loop is knotted across the first 
vortex loop. Comparing Eqs. ^ and ([^ suggests that 
the Chern-Simons contribution is nonzero only if the two 
vortex loops are linked. 

We verify the above proposition as follows. For two 
nonintersecting curves 71,72 (which act as mapping from 
manifold Si (circle) to 3-d Cartesian space (the 3-d 
Euclidean space in our problem), the linking number is 
given by 


M 


linking 


47r 
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r2) ■ (dri X dr2) 

|ri - i'2P 


( 10 ) 


Taking the inverse Fourier transform, with the vortex 
current given by Jy = qdr and endowing the above ex¬ 
pression with an overall coefficient k, the above expres¬ 
sion corresponds to vortex current-vortex current inter¬ 
action 


7^ 7 S 

( 6 ) 

In this case, we have decomposed the phase field into 
vortexful and regular parts (p = (pvY <pR defined vor¬ 
tex loop current Jy = {l/{2TT))e^^''d^d^(p\i. Following 
boson-vortex duality transformation m and working in 


Tiinking[Tv] = y Tv(fc)^ (-KTre^i^aA:") Jy(-fc) (11) 

in Fourier space. This is precisely equal to the topological 
part of Eq.Q. This implies that the Chern-Simons term 
leads to the linking between vortex-loops, resulting in 
nonzero linking number, precisely as implied by Eqs. ([^ 
and (|^. 
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III. ANYON PHYSICS 

We now show how the field theory and the results in the 
previous section lead to the anyonic fractional statistics of 
the vortex loops excitations of the 2-d XY model Eq. ® . 
Without the Chern-Simons term, the XY model Eq.® 
has topological excitations in the form of vortex loops 
which, in the absence of destructive interference due to 
the Berry phase, can proliferate. The vortex loops can be 
considered as worldlines of Bose particles and the prolif¬ 
eration of vortex loops corresponds to a condensation of 
the bosons. This vortex proliferation implies disordering 
of the dual {(p) field, which gives rise to gapped mag¬ 
netic excitations. In the presence of the Chern-Simons 
term, vortex loops become anyons and a number of vor¬ 
tex loops must collect together in order to form a boson, 
and then be able to condense. One then can have either 
one of the two scenarios: i)The anyonic vortex loops are 
able to form bosons and condense, giving rise to a state 
with gapped magnetic excitations as the analog of FQHE 
state, ii) The anyons fail to condense, for which a state 
with gapped magnetic excitations will not occur. 

To show this technically in a simple picture, consider 
a system of two vortex loops with vorticities qi and 52 
respectively knotted to each other once (Minking = !)■ 
The partition function of this system has a contribution 
of the form 


ii-KKqi q 2 +qi ll 2 ^^Coulomb 


( 12 ) 


Apart from the Coulomb interaction, the contribution to 
the partition function given in Eq. (121 resembles that 


of linked wordlines of anyon system with its fractional 
statistics [18j . The first term in the exponent 'i'!TKqiq 2 
equals twice the statistical angle 0 under the exchange 
of two anyons realized by vortex loops with vorticities qi 
and q 2 , giving 0 = 2'KKqiq2. In 3-d Euclidean space-time, 
there exists non-trivial braiding statistics between loops 
m- Our system thus manifests anyonic loop statistics in 
3-d. With qi — q 2 = i for elementary vortex loop, we ob¬ 
tain 0 = 27r/t and therefore need A/'anyon = Stt/O = 1/k 
anyons to get a boson before the anyons can condense. 
The occurrence of anyon condensation (and that of spin 
gapped state) thus depends crucially on the magnitude 
of the Chern-Simons coupling k. The Coulomb interac¬ 
tion then determines the critical value of the parameter 
that characterizes any possible condensation transition 
for this interacting system. We will show a scenario from 
a microscopic model where k takes small enough values, 
leading to large A/knyon and as a consequence prevents 
anyon condensation and the occurrence of a state with 
gapped magnetic excitations. 


IV. MICROSCOPIC REALIZATION 


spin field coupling. This topological term is obtained 
upon integrating out the fermions, as summarized in Eq. 
0 - The key ingredient is to have massive Dirac fermions 
coupled to the spin field from the spin sector for which 
we propose to consider the staggered tt flux state m and 
we will now describe the scheme to obtain it. 

According to Eq. ([^ , the bare hopping integral is mod¬ 
ified by the overlap of spin coherent states t —>■ t(f2r|Dr'). 
This spin coherent state overlap generally takes complex 
values which gives an exponential phase factor where the 
argument of the exponent, to be referred to as link spin 
field, contains a static background flux plus the fluctu¬ 
ating part, representing a dynamical spin field. In order 
to get Dirac spectrum, the static background flux has to 
take an appropriate configuration, which can be chosen 
to be staggered tt flux configuration [H]. In order to at¬ 
tain such staggered tt flux configuration, in turn, the spin 
sector should take a particular spin configuration accord¬ 
ingly. In the Appendix A, we first derive the expression 
for the link spin field from evaluating the spin coherent 
state overlap explicitly. Then, from that we determine 
the corresponding spin configuration to obtain the stag¬ 
gered TT flux state. We then propose a microscopic spin 
Hamiltonian that stabilizes the required spin configura¬ 
tion. We summarize the results as follows. 

Representing the spin as a classical vector S = 
S (sin 0 cos (j), sin 6 sin (p, cos 9) and evaluating the spin co¬ 
herent state overlap in Eq. ([^ give the following result for 
the link spin field [22] . 


= — tan 


(13) 

The flux $□ through a square plaquette is then given 
by $□ = X)rr'Gn®rr' where the link rr' is taken to be 
such that the four links traverse the square in counter¬ 
clockwise direction. The ’mean-field part’ of the link spin 
field needed to give tt-Aux state is obtained by assigning 
(pr their classical lowest-energy configuration values. 
It can be checked from Eq. (13), that a ±7r flux through 


square plaquette can be realized by planar spin config¬ 
uration where 0 = 7r/2 at all sites while the azimuthal 
angle takes values such that the four spins on a square 
rotate by a total of ±27r as illustrated in Fig. 0. The 
need for such planar spin configuration fully agrees with 
the intuition that the total solid angle $ swept by the 
four spins around the square should equal 27r for them to 
produce net flux $5 = 5'$ equals tt with S = 1/2. 

Similar planar spin configurations can be realized us¬ 
ing Ji — J 2 model with impurity [20]. Here we propose 
an alternative method. We find that such kind of spin 
configurations can be stabilized as the classical lowest- 
energy configuration of the following Hamiltonian, in¬ 
volving Dzyaloshinskii-Moriya (DM) interactions [2T]. 


In this work, we propose a way to generate a Chern- 
Simons term explicitly in a spin system from the fermion- 


77 = J^S.-S,+d5](5/)2+^D/^^.(S.xS,) (14) 

ij i ij 
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Eq.0 with V'(r) = {cjj{r),cl,{r),c%{r),c%{r))'^. Here 
Cp,i = I,-- - ,4 represents the Dirac fermion operator 
for the sublattice, while 7 ^, ^ — 0,1,2 are 4x4 ma¬ 
trices satisfying the Clifford algebra { 7 ^, 7 ;^} = and 
ijj = ipl'yQ [55]. We take the unit hopping integral t = 1 
and unit lattice spacing a = 1 as units of energy and 
length respectively. One important result is that the 
gauge charge is found to be 


FIG. 1. The planar spin configuration on square lattice with 
the labelling of the four sublattices. 


The easy-plane anisotropy {D > 0) favors in-plane Neel 
state. To retain global C/(l) symmetry, only the com¬ 
ponent of along z direction can be allowed to be 

nonzero. DM interaction usually occurs in lattice systems 
with low degree of symmetries. Here, we consider stag¬ 
gered DM interaction where 

while = 0 According to Moriya’s 

symmetry analysis |21] . this DM pattern is in principle 
allowed as long as the system has no bond-centered in¬ 
version symmetry but has mirror planes perpendicular to 
the bonds and passing through the bond centers. 

To find the lowest energy spin configuration, we re¬ 
write the Hamiltonian Eq.(14) in terms of angular vari¬ 
ables 0, (j) as before and minimize H with respect to these 
variables [55]. The lowest energy configuration of H in 
Eq. (141 is found to be a planar spin configuration with az¬ 
imuthal angles = (j)j — (j)i between nearest-neighbor 
spins given by 


= tan ^ 



(15) 


as shown in Fig.Q. Our calculation confirms that the 
planar spin configuration is indeed the lowest-energy con¬ 
figuration of H at least within a finite regime of phase 
diagram defined in parameter space {J—D — Dom space) 
and therefore requires no fine tuning. 

We then derive from Eq.([^ the low-energy effective 
theory of fermions around the Dirac points coupled to 
the spin quantum fluctuations around the lowest energy 
spin configuration. This is done by expanding the spin 
vector around its lowest energy orientation and expand¬ 
ing the fermion field operator using gradient expansion 
in real space and expansion around Dirac point in mo¬ 
mentum space. The details are given in Appendix B. The 
resulting scalar theory takes the form given in Eq. (|^ for 
the spin sector with (j) = (j)G = <)>i)/^unitceii 

(where (j)i is the (j) of the sublattice) representing 
the massless Goldstone mode. Here, as the other fields 
are gapped, gapped magnetic excitations imply spin gap. 
The Berry phase term gives a contribution exp(f27rS') 
to the partition function in imaginary time. For inte¬ 
ger S this equals unity corresponding to constructive in¬ 
terference and one can thus expect the presence of spin 
gapped state. The fermion-spin field theory is given by 


e 


1 

2^2 





(16) 


where a = e^/(47r) gives the dimensionless ’fine structure 
constant’ of the 17(1) gauge theory. This gauge charge is 
obtained directly as the coupling constant of the derived 
fermion-spin field theory. It is to be noted that in realistic 
situations, \D]:,m\ ^ J and as a result, a becomes a 
small parameter, as is normally the case in quantum field 
theory. It will be shown in the next section that this 
fact will play a key role in the mechanism of protection 
against spin gap that we propose in this work. Strictly 
speaking, the derivation of our fermion-spin field theory 
is justified only in the large spin S limit. While we have 
performed explicit calculation to obtain the quantitative 
result for the gauge charge e as above using S' = 1/2 [55] 
as example, our main result regarding novel mechanism 
for protection against spin gap that follows soon will still 
hold qualitatively and apply to general spin S. 

Now we show how to get the Dirac mass for the mas¬ 
sive Dirac fermion spectrum. We consider a perturba¬ 
tion on the ideal system in terms of dimerization of the 
strength of the hopping integral and staggerization of the 
sublattice chemical potential. Since we aim for a Chern- 
Simons effective theory, which intrinsically breaks dis¬ 
crete symmetries, the full system consisting of the spin 
and fermion sectors including the perturbation should 
break the symmetry under parity and time reversal times 
any lattice translation. Eq. Q also suggests that Dirac 
fermion masses from different fermion flavors must neces¬ 
sarily have the same sign in order to ensure nonzero net 
Chern-Simons term, as masses of opposite signs will give 
no net parity-breaking effect [53]. We find that this can 
be achieved by considering a combination of a columnar 
dimerization and a staggered sublattice chemical poten¬ 
tial with profiles shown in Fig.Q. It gives a Dirac mass 
Hamiltonian 


Hm = J (17) 

where V’i( 2 ): V’i( 2 ) represent the first (second) pair of sub¬ 
lattice l-sublattice-2 spinor ^’ 1 ( 2 ) = ('- 1 ( 3 )!‘' 2 ( 4 ))'^ 
agonalized basis, with 7 ° = r^, the Pauli matrix in 
this sublattice space. The resulting mass is found to be 
= V’^dimer + where rndimer = r? - 1 with ry > 1 
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FIG. 2. The columnar dimerized plus staggered chemical po¬ 
tential added to the staggered tt flux state. 


is the strength of dimerized bond with respect to nor¬ 
mal bond and rrig is the staggered sublattice chemical 
potential. The above perturbation is an example but in 
principle, any perturbations breaking the same discrete 
symmetries are expected to give rise to the same phe¬ 
nomena. 


V. DISCUSSION 


A Chern-Simons term is known to give rise to a mass 
to the gauge field in the Maxwell-Chern-Simons theory. 
What we have here is an ATT-Chern-Simons theory plus 
the ever important Berry phase term. Plateaus in XY 
model without a Chern-Simons term occur when the 
spin sector is gapped and vortex-loops proliferate cor¬ 
responding to condensation of bosons. In the presence 
of the Chern-Simons term, vortex loops are anyon world¬ 
lines as we noted from Eq. (12), with statistical phase 
0 = 2ktt = 12. If the vortex loops could ever con¬ 

dense, we would realize the analog of a fractional quan¬ 
tum Hall state (FQHE) with its chiral edge states. For 
our microscopic model with gauge charge e as given in Eq. 
(16), we obtain an upper bound 0 = 7r/16 from the limit 
\Ddm\I J oo which means that in order to obtain con¬ 
densation of anyons, we need at minimum ACanyon = 32 
vortex loops to form a boson first before they can ever 
condense and give rise to spin gapped state, taking the 
fact that for bosons, the statistical angle is 0 = 27r. Real¬ 
istic situations where \Dy)m\IJ “C 1 require much larger 
number of anyons for them to condense. 

Based on the analogy with the FQHE, condensation 
of anyons and thus Chern-Simons-induced spin-gapped 
state in doped antiferromagnets can occur only when the 
microscopic parameters give rise to a Chern-Simons term 
with coupling k such that 1 / = 2k = Q/7r = e^/2 = P/Q 
satisfies a continued fraction expansion condition, where 
iz is the filling factor of typical values l/v < 9 for Laugh- 
lin’s states [24]. With a minimum of 32 vortex loops 
needed to form a boson (equivalent to 1/u 16 for the 

realistic case \Ddm\/J ^ 1) and the strict condition 
on the Chern-Simons coupling, it is thus in general very 
unlikely to form such vortex loops condensate and the 
associated spin gapped state. We can therefore conclude 
that the spin system Eq.( 14) is protected from being in a 


spin gapped state, due to the Chern-Simons term induced 
by the fermion-spin field coupling in doped antiferromag¬ 
nets with a massive Dirac fermion spectrum. In contrast, 
the conventional XY model with the Berry phase term 
would mandate the occurrence of a spin gapped state. 
On the other hand, if the net Chern-Simons term van¬ 
ishes, then the protection effect is inactive and one can 
get back a conventional spin gapped state that occurs 
for integer spin S. It is to be noted that this result on 
the novel mechanism for protection against spin gap is 
valid for general spin S because the applicability of the 
mechanism is determined more by the value of coupling 
constants of the spin model (J and Ddm) rather than 
the spin S itself. 

It would be interesting to find other microscopic spin 
models which can generate Chern-Simons term via sim¬ 
ilar fermion-spin field coupling as we proposed here and 
yet are able to induce anyonic vortex loop condensation, 
giving rise to a spin gapped state, and the analog of the 
FQHE with chiral edge states in spin systems. Our choice 
for square lattice is because a discrete Chern-Simons 
gauge theory with precisely the same physics as contin¬ 
uum one can be constructed consistently on this lattice 
pS] as well as on kagome lattice [15] , due to the presence 
of one-to-one face-vertex correspondence on these lattices 
[27] . We would like to propose a study on compounds 
La 2 Cu 04 [28] and LaMnOa [29] which are effectively 2-d 
square lattice antiferromagnets as the candidate materi¬ 
als to test our theory. Beyond these particular materials, 
we claim that this novel scenario of spin gap protection 
is general enough to be expected anytime a Chern-like 
charge insulator is at play as the high energy sector. The 
Kondo lattice model in the triangular lattice may be an¬ 
other possible laboratory [30] . 
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Appendix A: The Classical Lowest-Energy Spin 
Configuration 

A key requirement in our proposal to generate Chern- 
Simons term in a spin system is to have Dirac spectrum 
for the fermions that are to be doped into the system. 
In order to get a Dirac spectrum on the square lat¬ 
tice, we need a tt-Aux state of the type as first proposed 
by Affleck-Marston and G. Kotliar m- The flux per 
square plaquette is a staggered tt and —tt configuration, 
where the flux alternates from tt to —tt between adjacent 
squares. Here we will show that such staggered tt flux 
state can be realized by a planar spin configuration illus- 















6 


2 _ 















- i 


-- 


--i. 


FIG. 3. The (commensurate) planar spin configuration on a 
square plaquette of the square lattice giving rise to staggered 
TT flux state of Affleck-Marston ansatz. 


trated in Fig.(|^. The derivation to arrive at such spin 
configuration is as follows. 

The spin sector couples to the fermion by providing 
the spin field in the form of Goldstone modes in addition 
to providing the background staggered rr-flux state. This 
is well described by the following tight-binding Hamilto¬ 
nian. 


i? = —t -b/l.c. (Al) 

r,r' 

In spherical coordinate, the classical spin vector is writ¬ 
ten as S = S {sin 6 cos (f), sin 9 sin (p, cos 9). The spin co¬ 
herent state for general spin S is given by 




mz — — S 


(cos )‘S'+”j.>!(sin '■*2 

y/{s + m^)l{S - m^)l 


-\S, m^) 


(A2) 


r 


which for 5" = 1/2 allows us to write 


|r2r) = i 


_ _. dj. 

e" 2 cos —, e 2 sin — 
2 2 


(A3) 


where is a pure gauge function which we can set to 
constant function = 60 in the simplest case. The co¬ 
herent state overlap for electron where S = 1/2 is given 
by 


(f2r|Hr') = e2 


— „5df'r'+dr'-^'r-0r) 


dr' 

cos — COS — 
2 2 


pi(0r-0r') 


dr dr' 

sm — sm — 
2 2 


(A4) 


This spin coherent state overlap gives rise to a spin field can be obtained as the lowest energy configuration of the 

as follows. Taking the gauge function to be 6 r = —pr, following Hamiltonian, 


, sin(d)r' — Pr) sin % sin -p- 

= -tan-^ --AT- — -^- h -§T 

cos ^ cos + cos((/r' — Pr) sin ^ sin 

(A5) 

We will consider the Dirac spectrum obtained from the 
staggered tt flux state configuration of Affleck-Marston. 
It can be checked from Eq. (A5) that for dr = 7 r /2 at all 
sites r,r', we will always get ^Jrr'en “^r' = as the 
link spin field is then given by Orr' = — Af/rr '/2 where 
A(/rr' = Pr' ~ Pr- This gives rise to the needed tt flux per 
square plaquette. This remarkable result will give rise 
to the planar spin configuration shown in Fig.(|^ for an 
appropriate choice of Hamiltonian. 

We will prove here that this planar spin configuration 


= jESrS,+DE('5fr+El^i^M-(SzxS,) (A 6 ) 

ij i ij 

consisting of a Heisenberg antiferromagnetic coupling 
(J > 0), easy-plane anisotropy {D > 0), and a 
Dzyaloshinskii-Moriya terms. To retain global U{1) sym¬ 
metry, only one {x, y or z) component of normal 

to easy plane can be allowed to be nonzero; for concrete¬ 
ness we choose 2 ; component. It will be shown here that 
in order to get the planar spin configuration Fig. (§ , we 
must take > 0 on cc links and < 0 on y 

links (of equal magnitude but opposite in sign) or vice 
versa. 

In spherical coordinate, the full Hamiltonian is given 

by 


H = JS^ [sin 9i sin 9j cos{pj — pi) + cos 9i cos 9j] -b DS^ E ^os^ d^ 


'ij 


E sind^ sin9j sin((/j - pi) (A7) 


To find the lowest-energy configuration, we take dH/dpi 


0, dH/d9j 


0 and verify whether 






















7 


d^H/d(l)id(l>j > 0,d^H/de,d9j > 0,d^H/d(j)id9j = 0 . 
We find the lowest energy configuration solution where 
the nearest-neighbor spins have polar angles as 9^ = 7r/2 
at all sites i and azimuthal angle Atjjij = (j)j — (j)i between 
nearest-neighbor spins given by 


A(j)ij = tan 


-1 / ^DM 


(A8) 


whereas next-nearest neighbor spins across the diagonal 
of square have A(j)ij = 4>j — 4>i = as shown in Fig.([^. 
We verify that the conditions on the second derivatives 
are indeed satisfied. The above calculation guarantees 
that such planar spin state indeed exists within a finite 
regime of the phase diagram defined in the parameter 
space {J — D — Djjm space) and therefore requires no 
fine tuning. 

Now we can describe the precise resulting lowest en¬ 
ergy spin configuration for a particular pattern of the 
Dzyaloshinskii-Moriya interaction coupling of rele¬ 

vance to our purpose. For case 1), we take = TDdm 
to be uniform equal to a constant vector at all links, and 
without loss of generality we take it to point to the z 
direction for simplicity. It can be checked from Eq. (A8) 


that unless {TAdmI/J = 1 giving = 7r/4, 57r/4, one 
cannot get tt flux on a square plaquette for realistic case 
where |D£)m|/< 7 ^ 1- The four spins just cannot be ar¬ 
ranged to rotate in total by 27r and at the same satisfy 
Eq. (A8| from link to link on a square. What this means 


is that we have to consider a plaquatte involving further 
distanced spins to get tt flux. Thus in general we will 
get an incommensurate state for case 1). For case 2), we 
consider lL>z)Mi,i±x = (0; 0) on each ’column’ and al¬ 

ternating lL)z)Mi,i±y = ( 0 ,±(i 2 ) 0 ) on each ’row’ of square 
lattice, where |di| = \d 2 \ and so di = ±d 2 - It can be 
easily checked using Eq. (A8) that one can arrange the 


four spins to rotate by 27r on a square but when we tile 
up the pattern on a lattice, it will take a relatively large 
cluster of squares with their spins to form unit cell from 
which the full lattice spin configuration can be formed. 
We however get a commensurate staggered tt flux state. 
This DM pattern occurs in real materials, e.g. La 2 Cu 04 
pS] but such large unit cell complicates the the deriva¬ 
tion of Dirac fermion field theory as the spinor and ma¬ 
trix sizes grow with unit cell size. We then consider the 
model case 3) where we have 

while = 0. It can be checked using 

Eq. (A8) that the four spins on a square will take con¬ 
figuration as shown in Eig. (|3) and this can be nicely 
repeated over the whole square attice and gives what we 
call (commensurate) planar spin configuration with the 
resulting commensurate staggered tt flux state with just 
four sublattices (four sites per unit cell). The needed DM 
pattern is in principle allowed by symmetry according to 
Moriya’s original consideration. We take this latter case 
as example but our final results apply to the real mate¬ 
rials mentioned above. 

We label the four sublattices as shown in Eig.Q. 



FIG. 4. Labeling of the four sublattices is as shown in the 
square. The labeling of the four Dirac points in Brillouin 
zone also follows the same numbering. 


It is to be noted that in such planar spin con¬ 
figuration, each sublattice has its own 9 and </>. 
We therefore have four sets of {(j)^''\9^‘^'>) fields: 

From these, 

we can define four orthogonal fields; 


4 >g = ^ { 4>i -\- 4>2 + 4>^ + 4>i )! 0 ti = ^ {4>i + 4>2 — — (i>i) 


4>t2 = ^ {4>1 — 4'2 + (t>3 — M ^ {(t>l — 4’2 — 4>3 + (j)4) 

(A9) 

~ ^ (^1 + 92 + 9^ + 94 ) ,9ti = - {9i + 92 — 93 — 94 ) , 


dt2 = ^{dl — 92 + 93 — 94 ) , 0(3 = - (01 — 02 — 03 + 94 ) 

(AlO) 

and the inverse mappings 

</’! = {4>G + (t)tl + 0(2 + 4>t3) , (t>2 = {(t>G + 4>tl — (t>t2 — 4>t3) , 


4>3 = (</>G — ^il + </>t2 — ^(> 43 ) ) </>4 = ((/'G — — 4’t2 + (j)t3) 

(All) 

= (^G + dtl + 0(2 + 9t3) , 02 = {9g + 9ti — 0(2 — 9t3) , 


^3 — (^G ~ ^(1 + ^(2 — ^ 43 ) ) ^4 — [da — 9ti — 0(2 + ^ 43 ) 

(A12) 

The effective low energy theory of the spin sector from 
the Hamiltonian Eq. (A61 will give rise to a gapless La- 


grangian for (pc and gapped ones for 0((, i = 1, 2, 3, so 
that these latter fields can all be integrated out in the low 
energy physics m- We therefore have only the (jra re¬ 
maining in the low energy physics. This (pc is nothing but 
the Goldstone modes of the rotationl symmetry-breaking 
classical lowest-energy configuration that we have ob¬ 
tained in the beginning: the planar spin configuration. 
The effective low energy description of the spin sector 
obtained upon integrating out all massive fields take the 
form 











demonstrate in this section the way to generate massive 
Dirac fermion spectrum. 

To derive the low-energy effective theory of fermions 
coupled to the spin field, we start from Eq. (All. 
Fermion in the staggered tt-Aux state of Affleck-Marston 
has Dirac spectrum Ek = ±2xiy^coPXT+cos^"^ with 
four Dirac points at k = (±7r/2, ±7r/2) as shown in 
Fig .([^ with their numbering. In the 4-sublattice descrip¬ 
tion of Affleck-Marston tt-Aux state, the new Arst Bril- 
louin zone is deffned by the square with kx-, ky = ±7r/2. 


FIG. 5. The Brillouin zone of the staggered tt Aux state in We Arst expand the fermion held operator around these 
the 4-sublattice desciption and the Dirac points (solid circles). four Dirac points 


C = ^d^<j)Gdf,(j3G + i^dr(l)G (A 13 ) 

and it is to be noted that the Berry phase term is ob¬ 
tained from summing over the Berry phases of the four 
sublattices which Anally involves only the diagonal held 
(Goldstone mode) with overall prefactor S. Since 
upon integration over r Berry phase gives phase factor 
exp(i27rS'), it equals unity for integer spin S and thus has 
no effect in such case. To simplify notation, in the rest 
of this calculation, (j),9 represent 0 g,^g- 


Appendix B: Derivation of Effective Fermion-Spin 
Field Action of Insulating Antiferromagnet with 
Massive Dirac Spectrum 

We have shown in the previous section the microscopic 
spin model to obtain (gapless) Dirac spectrum for the 
fermions doped into the system. We have to verify that 
the resulting low-energy effective theory is indeed that 
of Dirac fermion. Furthermore, in order to generate the 
Chern-Simons term properly, the Dirac fermions must 
be massive rather than massless. We will therefore also 



= (Bl) 


where the slow-fermion field operator CDa(j) is given by 


CDa(r) 



d^k 


Cke 


z(k—koa)- 


r 


d^q 


'|qa|<9< 


(27r) 


2 ^^na+q.a^ 


lC\a 'T 


(B 2 ) 


with Qa = k — and \i.Da are the momentum of the 
Dirac points \^Da = (±7’‘/2, ±7r/2) and qc is some 
appropriate UV cutoff. 

We then perform gradient expansion and in doing so, 
we expand the spin coherent state parameterized in terms 
of these two angles around the respective classical lowest- 
energy conAguration value 6q and at each site (sublat¬ 
tice). Eventually, only the Goldstone mode (pG survives 
and in the following, sublattice index is omitted and in 
the Anal results (j),9 represent Pg^Og- 


|Or( 0 r,<^r)) = \^M. €)) + . €)) + . €)) + 0{69\ 5cp\ 595^) (B 3 ) 

I {9r'Ar')) = \ {0 °', )) + Off \ Dr- {9°,, P°,)) + 5cj>x' 0^ \ Dr- ( 0 °,, P°,)) + 0{69 ^, 5 ^^, S9SP) (B 4 ) 


■ '^^Da + 0{Ar^) (B5) where Ar = r' — r, we obtain from Eq. (A1) 


r.r' a,b—l 


;(kDb-r-koa-r')r„t rryoOlOO 


CDai'^)i^r\^r')cDbir) - • (-t(Dr|Dr-)V - eA)cDb{r)] + h.c. 


l^r+Ar)^-C>a(l‘) ' i ^(^rl^r+Ar)^ €.A.)cjp}a{^)] h.C. (^6) 

r a—1 


.t 
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4 4 


= -^EE 


'y ' g-ikna-Arf^mt 


- c'"I,„*Ar • i-^{n^M+Ar)^ - eA)c£,(r)] + h.c. (B7) 


r a—1 m,n—l 


r 


where the superscript m,n = 1, • • • ,4 represents sublat¬ 
tice index of the spin sector and Da = 1, • • • ,4 the Dirac 
point index. In what remains, we take the unit hopping 


integral t = 1.0 and the unit lattice spacing a = 1.0 as 
the units of energy and length respectively. The spin field 
is given by 




(B8) 


where in the continuum limit we do gradient expansion 
by defining V0 = 69ri — 6dr,V4' = S(j)r' — Scj)r. We have 
included the gauge charge e couplin g th e spin field to 
the fermions. The expression Eq. (B8) suggests that 
this spin field is a Berry connection-type of pseudo-gauge 
field. The actual contribution to the spin field A in Eq. 
(B8| only comes from the (j) field part, corresponding to 
Goldstone mode representing symmetric combination of 
the sublattice phase fields (j)G = {(j)! + (/'2 + ^3 + </>4)/4, 
whereas the terms containing 59^, S9r' give rise to massive 
momentum operator and correspond to the massive 9 
field, which can therefore be integrated out in the low 
energy physics. 

Eor S = 1/2, we get the result for the spin field 


eA(r) = — ^ sin^o sin V(/(r) (B9) 

where 9o = 7r/2. The above result can be derived in a 
more direct and straightforward way by evaluating the 
coherent state overlap and matching the complex phase 
factor exponent with spin field. 


(Dr|Dr-) - e“--' (BIO) 

and as one takes the continuum limit, one has arr' —t 
A(r). We obtain 


(Dr|Dr') 






J./ -4>r-br ^ 



where = (pj./ — and 


= sin 00 cos 


A())rr' 

2 


(B12) 


where 0o = 7r/2 for our planar spin configuration. We 
can immediately identify 


eA(r) = — - sin 00 sin 




,r-|-Ar 


V</(r) 


(B13) 


in complete agreement with Eq. (B9), while the last term 
in Eq. (Bll) is dependent on the momentum operator 


Hr = — ^[<50rSin0r -I- ^50r COS0r] (B14) 

which is massive and therefore in the low energy physics 
can be integrated out to give simply a renormalization 
correction to the coefficients (couplings) in the effective 
fermion-spin field action without changing the physics. It 




1 ^ ^ / 

-]—sin —sinSoV^+i- 


/ 

- tan —(rTr+rij./) 


(Bll) 


is to be noted that the expression for spin field Eq.(B13) 


is proportional to the gradient of the phase field, as we 
expected. _ 


Following our definition = d^(j) and using Eqs. (A8) 


coupling 


and (B13), the above result gives as the effective gauge 


e 





1 

2 


(B15) 


We note that if we set the Dzyaloshinskii-Moriya inter¬ 
action to zero Ddm 0, the whole picture breaks down 
and the gauge charge equals zero. So, DM interaction is 
crucially needed here. 

As can be seen from Fig.([^, the Dirac points k^a = 
(±7r/2, ±7r/2) are located precisely at the corners of the 
fst Bi-iPouin zone. The implication of this is that the 
four Dirac points are now equivalent to each other and 
it suffices to consider only one Dirac point, e.g. the 
koi = ('^/2,7r/2) Dirac point. In the low energy limit. 
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we define a Dirac 4-spinor associated with this represen¬ 
tative Dirac point '0(r) = (c}j(r), c|,(r), c|j(r), c|j(r)) . 
Here, c} 5 (r), i = 1, • • • ,4 represents the Dirac fermion an¬ 
nihilation operator for the sublattice of the spin sec¬ 
tor in the planar configuration Fig. The low-energy 
Hamiltonian is 


H = 


- eA^)tp 


(B16) 


where /i = t, a:, ?/ = 0 , 1 , 2 and ip = ip^'^j^. 

The next task is to find the 7 matrices. We obtain the 
following result. 


7 = *7 


0 7^2 0 0 \ 

7fi 0 0 0 

0 0 0 7^4 

0 0 7f3 0 / 


(B17) 


y • t 


000 7^4 
0 0 7 I 4 0 

0 7f2 0 0 

7I1 0 0 0 


(B18) 


where the matrix elements are given as follows. 




7^4 = ite = ( 7 I 4 ) 


= it 


with 


7 = 


10 0 0 
0-100 
0 0 10 

0 0 0 -1 


= (7*)- 


(B22) 


We aim this Dirac fermion system to produce a Chern- 
Simons term. For that we need the spectrum to be mas¬ 
sive. We need to make sure that the Dirac masses from 
the first two sublattices and the second two sublattices 
have the same sign, otherwise the net Chern-Simons term 
vanishes. Such mass term can be obtained in several 
different ways. First, using staggered sublattice chemi¬ 
cal potential which represents some internal sublattice or 
microscopic degree of freedom which gives rise to effec¬ 
tive chemical potential alternating in sign between the 
four sublattices of the square lattice. In this case the 
mass term from staggered chemical potential takes the 
form Hm = / [diag(ws, —mg, TOg, —TOs)]V’ where 

[diag(- • •)] represents diagonal matrix with the diagonal 
elements as given. 

The second method is to consider dimerization of the 
strength of the hopping integral t in the tight-binding 
Hamiltonian Eq.(Al|. Dimerization has been known to 


produce a mass for the Dirac fermions. We consider one 
of the simplest dimerization patterns as shown in the 
Fig. (§. The strength of the dimerized bonds is larger 
than the undimerized bonds. The link spin field and 
the flux per plaquatte remain the same as those of the 
original Affleck-Marston tt-Aux state. This dimerization 
does open up gap on the Dirac points in the 4x4 form 
of the Dirac theory which for the dimerization pattern 
shown in Fig. (§ takes the form 


7 I 3 = = (7|2)^ 

7-4 = = ( 74 - 3 )^ (B19) 


/ 0 mo 0 0 \ 

o/dimer _ t ^0 ^ ^ ^ I 

- 1 0 0 0 mg 

\ 0 0 mg* 0 / 


(B23) 


These matrix elements need to be ’normalized’ in such a 
way that each of them is of unit magnitude. We have to 
find the appropriate 7 ° = 7 *. It needs to be checked that 
the 7 matrices satisfy Dirac algebra { 7 '^, 7 '^} = 25^'^. 

Taking t = 1, we obtain the following gamma matrices 
valid for a spin system in planar configuration on the 
square lattice. 


where mo = —(77 — l)exp(—i7r/4 — ikux),mQ = —(77 — 
1 ) exp(—i 7 r /4 -|- ikox) with 77 > 1 is the relative strength 
of dimerized bond with respect to undimerized bond of 
unit strength. Using kox = 7r/2 from the location of the 
Dirac point, we obtain mg = —mo = ( 77 —1) exp(—377r/4). 
So, the Dirac mass Hamiltonian induced by dimerization 
Fig. ([§) can be written as 


7 =77 



1 

0 

0 \ 

(0 

7 

0 


1 1 

0 

0 

0 1 

1 _ I -* 

0 

0 

0 

0 

0 

0 


1 “ 0 

0 

0 


Vo 

0 

-1 

0 ) 

V 0 

0 

7 

0 / 


(B20) 


7 


y 


= i"f 


t 


/O 0 0 
0 0 1 
0 1 0 
V 1 0 0 



0 0 0 7 \ 

0 0 -7 0 I 

0 7 0 0 

-7 0 0 0/ 


(B21) 


Hm = 


d^xiputr^i^ = 




(^-1), 
72 


T'^)a^ip 

(B24) 


where r is the Pauli matrix defined within each of the 
1 — 2 and 3 — 4 sublattice pairs while a is the Pauli matrix 
connecting the two pairs. In the rest of this derivation, 
we always imply 


Hm = 


cPxip'H.M'ip 


(fxip^WM'tp 


(B25) 
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where Hm = ( 7 *) ^T-Lm is the Hermitian mass matrix 
derived directly from the tight-binding Hamiltonian. 


TT 

^ i 

— TT 

IT 

— TT 

TT 

— TT 

TT 

— TT 

TT 


FIG. 6. The columnar dimerized tt flux state of Affleck- 
Marston ansatz. 


Evaluating the mass Hamiltonian matrix Eq. (B23) at 
the representative Dirac point, we obtain the following 
result for mass matrix. 


■T/dimer 

tt-M — 7 


/ 0 


V 


(’?-!) 

72 


( 1 -i) 


0 

0 


(^- 1 ) 

72 


(1 + *) 


0 

0 

0 


0 

0 

0 


72 


(1-z) 


0 \ 
0 

0 y 


(B26) 


The above result however suggests that the masses are 
opposite in sign between the two ( 2 -sublattice) pairs of 
the 4-sublattice theory (1-2 sublattice pair versus 3-4 sub¬ 
lattice pair). The net Chern-Simons term will be zero 
in this case. To remedy this, we combine columnar or 


staggered dimerization with staggered sublattice poten¬ 
tial where the staggerization takes 2 -sublattice structure; 
the chemical potentials of sublattices 1 and 3 are equal 
but opposite to the chemical potentials of sublattices 2 
and 4. The resulting mass matrix is 


/ 


/T/dimer+stagg.chem.pot 


= 7 


V 


ms 


(^- 1 ) 

72 


(1 


0 

0 


i) 


-^(1 + *) 0 0 \ 

—ms 0 0 

0 ms + 0 

^ -'ms ) 


(B27) 


where ms is the mass from staggered sublattice chemical 
potential. It can readily be seen that the above mass 
matrix will give rise to nonzero net Chern-Simons term 
since the mass block-matrix of the 1 — 2 sublattice pair 
is not opposite of that of the 3 — 4 sublattice pair as the 
staggered sublattice potential-induced mass term matrix 
has precisely the same form and sign in both sectors. If 
upon diagonalization the masses are of the same sign in 
both sector, which is indeed the case at least for certain 
range of parameters rrig, the the net Chern-Simons term 
is nonzero. Performing this 4x4 matrix diagonalization 
of the corresponding Mm, we find that the Dirac mass 
Hamiltonian can be written as 


where '*/'i( 2 )) V'i( 2 ) represent the first (second) pair of sub¬ 
lattice l-sublattice -2 spinor ijj'i = ( 0 ^, 02 )^ in diagonal¬ 
ized basis, with 7 ° = and mass = \/+ m^ 
where mdimer = rj — 1 for both pairs, indicating that 
the two pairs have identical mass and thus nonzero net 
Chern-Simons term eventually. This result is nicely con¬ 
sistent with symmetry considerations because the combi¬ 
nation of columnar or staggered dimerization plus stag¬ 
gered chemical potential completely breaks parity. It 
can also be checked that this dimerization plus staggered 
chemical potential breaks the invariance under time re¬ 
versal times any lattice translation and this guarantees 
nonzero net Chern-Simons term. 


-Hm = 


-I- 'ip'2m,p'ip2] 


(B28) 
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